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ABSTRACT 

Recent analysis revealed that supergranules (convection cells seen at the Sun's surface) are 
advected by the zonal flows at depths equal to the widths of the cells themselves. Here we probe 
the structure of the meridional circulation by cross-correlating maps of the Doppler velocity signal 
using a series of successively longer time lags between maps. We find that the poleward meridional 
flow decreases in amplitude with time lag and reverses direction to become an equatorward return 
flow at time lags > 24 hours. These cross-correlation results are dominated by larger and deeper 
cells at longer time lags. (The smaller cells have shorter lifetimes and do not contribute to the 
correlated signal at longer time lags.) We determine the characteristic cell size associated with 
each time lag by comparing the equatorial zonal flows measured at different time lags with the 
zonal flows associated with different cell sizes from a Fourier analysis. This association gives 
a characteristic cell size of ~ 50 Mm at a 24 hour time lag. This indicates that the poleward 
meridional flow returns equatorward at depths > 50 Mm - just below the base of the surface 
shear layer. A substantial and highly significant equatorward flow (4.6 ± 0.4 m s^^) is found 
at a time lag of 28 hours corresponding to a depth of ~ 70 Mm. This represents one of the 
first positive detections of the Sun's meridional return flow and illustrates the power of using 
supergranules to probe the Sun's internal dynamics. 

Subject headings: convection, Sun: rotation 



1. INTRODUCTION 

Supergranules are cellular flows observed in 
the Sun's photosphere as a s urface filling pattern 



of la rgely horizontal flows (jRieutord fc Rincon 
20101) . Typical supergranules have diameters of 
~ 30 Mm and lifetimes of ~ 1 day - much larger 
and longer lived than granules (diameters of ^ 1 
Mm and lifetimes of 5 minutes) . 

Since their discovery by iHart ( 1954) and 
their initial characterization by iLeighton et ah 



(|1962I ). supergranules have been recognized as 
the principal organizer of the Sun's chromo- 
spheric/magnetic network and as key to the 
transport of magnetic flux in the photosphere 
(jLeighton 1 ll964t iDeVore et al. lIlQSl . 

The rotation rate of the supergranules has been 
measured by a n umber o f investigators with in- 
teresting results. iDuvall I (|l980l) cross-correlated 



the Doppler velocity pattern from equatorial spec- 
tral scans obtained over several days and found 
that the pattern rotates about 3% faster than the 
photospheric plasma. Faster rates are found for 
the 24-hr time lags from day-to-day than for the 
8-hr time lags f rom the beginn i ng to end of an 
observing day. iBeck fc Schou ] (I2OOOI ) measured 
the rotation of the Doppler velocity pattern us- 
ing a Fourier technique and found that the larger 
cells rotate more rapidly than the smaller cells 
but with apparent rotation rates that exceed the 
peak internal rotation rate. This "superrotation" 
was attributed to line-of-sight pr ojection effects on 
the D oppler velo city pa t tern b y iHathawav et al. 



(|2006l ). In fact, ISchou I (|2003l ) had earlier found 
rotation consistent with the internal plasma flow 
proflle when this line-of-sight projection was re- 
moved. 

These observations are consistent with the orig- 
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inal conclusion bv lDuvalll ()l980l) that larger cells 
dominate the longer time lags and that these 
larger cells are more deeply anchored in a surface 
shear layer in which the rotation rate increases 
with depth. While it should be noted that some 
characteristics of the motions of the supergran- 
ules have been attr i buted to wav e-like properties 
(jGizon et al. II2003L ISchou II2003I ). it is clear that 
the cellular structures are advected largely by the 
zonal flows in the surface shear layer. 

The ex istence of this sh e ar la yer was sug- 
gested by iFoukal &: Jokipii ( 19751 ) as a conse- 
quence of the conservation of angular momentum 
by convective elements moving inward and out- 
ward in the near surfac e layers. Helioseism ology 



(jThompson et al. I Il99d : ISchou et al.1ll998l) 



indicates that the shear layer extends to a depth 
of ~ 50 Mm. 

Recentlv iHathawav ( 2012 ) found a one-to-one 
correspondence between the rotation rate of su- 
pergranules with increasing size and the rotation 
rate in this surface shear layer with increasing 
depth. Results, using a Fourier technique like that 
of , Beck fc Schou I (|200fl[ ) with the line-of-si ght pro- 



j ection effects removed using the method of lSchou 
(|2003| ). indicate that supergranules are advected 



by the flows in the shear layer at depths equal to 
their widths. We expect the advection of the su- 
pergranules to be dominated by flows near their 
bases since the plasma density increases rapidly 
with depth. However, this does imply that the 
cells must extend even deeper than this anchoring 
level. 

Note that while cells that are as deep as they 
are wide are typ ical in numerical simulations of 



solar convection (jStein et al. ll201lD . results from 



local helioseismology have suggested flatt ened su- 
pergra nules - but with conflicting results. iDuvall 



(|1998[ ) found that the outflows in the supergran- 
ules reversed at a depth of 5 Mm and then dis- 
appea red by 8 Mm depth. IZhao fc Kosovichev 



(|2003l ) found reversal at a depth of 8 Mm a nd dis- 
appearance at 15 Mm. ISvanda et al. 1 (l2009l) found 
disappearance at an even greater depth - 25 Mm. 

Here we measure the meridional motion of the 
pattern of supergranules by analyzing the sarn e 
dat a used inlBeck fc SchmTI (|2000[) . [Schou I (|2003f ) . 
and lHathawavl (|2012l) . 

The Sun's meridional flow is extremely weak 



when compared to the other photospheric flows. 
Both the strength and the direction of the 
meridional flow were uncertain until data and 
data an a.lysis techniques we re sufficiently im- 



proved dKomm et al. I 1993 ; Hathawav et al. 



1996l :lGiles et al. lll997l:rUlrich II2OIOI) . The merid 



ional flow is now known to be poleward in each 
hemisphere with a peak velocity of about 10-20 m 
s^^ while the zonal flows (differential rotation) are 
an order of magnitude stronger and the convective 
flows (granules and supergranules) are yet another 
order of magnitude stronger. 

The Sun's meridional flow plays a critical role 



in th e transport of magn etic flux (jDeVore et al 



1984 IWang et al. 1 l2009t ) and in some dynamo 



models for the Sun's magnetic cycl e (jPikpati &: Charbonneau 
19991 : iNandv fc Choudhuri 1 1200^. Although the 
nonaxisymmetric convective motions transport 
magnetic elements at much higher speeds to the 
cell boundaries, those motions are in random di- 
rections and give diffusion away from flux concen- 
trations. The net poleward transport has a larger 
contribution from the slow, but direct, meridional 
flow. This poleward transport is directly respon- 
sible for the reversal of the Sun's polar magnetic 
fields every solar cycle and for producing the polar 
fields that determine the strength of the following 



solar cycle (jSchatten et al. 
200,5h . 



19781 : iSvalgaard et al. 



Fourier techniques (which do give explicit in- 
formation on the cell sizes) are, unfortunately, not 
available for measuring the meridional motion of 
the supergranule Doppler velocity pattern. Here 
we use cross-correlations with time lags that vary 
from 1 hour to 28 hours. The displacement of 
the peak in the cross-correlation gives both zonal 
and meridional velocities which vary systemati- 
cally with time lag. We determine the charac- 
teristic cell size for each time lag by comparing 
the zonal velocities with those obtained with the 
Fourier technique ( Hathawav 20121) . This then 



gives the meridional flow at a series of anchoring 
depths corresponding to the series of time lags. 

2. DATA PREPARATION 

The data consist of 1024^ pixel images of the 
line-of-sight velocity determined from the Doppler 
shift of a spectral line due to the trace element 
nickel in the solar atmosphere by the Michelson 
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Doppler Imager ( Scherrer et al. 19951 MDI) on 
the ES A/NASA Solar and Heliospheric Observa- 
tory (SOHO). 

The images are acquired at a 1 minute cadence 
and cover the full visible disk of the Sun. We av- 
erage the data over 31 minutes with a Gaussian 
weighting function which filters out variations on 
time scales less than about 16 minutes, and sam- 
ple that data at 15 minute intervals. We then 
map these temporally filtered images onto a 1024^ 
grid in heliographic latitude from pole to pole and 
in longitude ±90° from the central meridian (Fig- 
ure 1). This mapping accounts for the position 
angle of the Sun's rotation axis relative to the 
imaging CCD and the tilt angle of the Sun's ro- 
tation axis toward or away from the spacecraft. 
Both of these angles include modifications in line 
with the most recent determinations of the orien- 



tation of the Sun's rotation axis (Beck fc Giles 



20051: iHathawav fc Rightmire II2010I) . We analyze 
data obtained during two 60 day periods of con- 
tinuous coverage: one in 1996 from May 24 to July 
22; the other in 1997 from April 14 to June 17. 

Each filtered ir nage is anal y zed u sing the meth- 
ods described bv iHathawav I ( 1992[ ) to determine 
and remove the stationary signals due to differen- 
tial rotation, convective blue shift, and the merid- 
ional flow - as well as the time varying signal due 
to the motion of the spacecraft. Following this, an 
average image from all of the images in each 60 day 
period was subtracted from each image to remove 
any instrumental artifacts. This leaves behind the 
non-axisymmetric and time-varying signal due to 
the convective flows. 

3. CROSS-CORRELATION ANALYSIS 

We determine the zonal and meridional veloc- 
ity of the supergranule pattern by cross-correlating 
strips of data 11 pixels high in latitude by 601 pix- 
els long in longitude. The strips are centered on 
each of the 860 latitude positions from 75° south 
to 75° north of the equator with overlap from adja- 
cent strips. These data strips are cross-correlated 
with corresponding strips from images obtained at 
time lags of 1, 2, 4, 8, 16, 24, and 28 hours. 

This process is done both forward and back- 
ward in time. A data strip centered on the cen- 
tral meridian is correlated with strips offset to the 
left at later times and then data centered on the 




Fig. 1. — Heliographic map detail of the line- 
of-sight (Doppler) velocity from SOHO/MDI. The 
map detail extends 90° in longitude from the cen- 
tral meridian on the left and about 35° in lati- 
tude from the equator (the thick horizontal line). 
The mottled pattern is the Doppler signal (blue 
for blue shifts and red for red shifts) due to the 
supergranules. 

central meridian at later times is correlated with 
strips offset to the right in the initial data map. 
This minimizes any systematic errors. 

The shift of the Doppler pattern is determined 
to the nearest pixel by finding the displacement 
of the strip in the later Doppler map that gives 
the maximum correlation. This pixel-wise shift is 
then refined by fitting parabolas through the cor- 
relation amplitudes at the displacement that gives 
the maximum correlation and those at plus and 
minus one pixel in each direction. The maxima 
of these parabolas then give the Doppler pattern 
shift to a fraction of a pixel. The shifts in lon- 
gitude give the apparent zonal velocity while the 
shifts in latitude give the meridional velocity. The 
60 days of data in each data set give some 1400 
independent measurements of these velocities at 
each of the 860 latitude positions for each time 
lag. 

The meridional flow measurements for time lags 
from 1 to 24 hours are shown in Figure 2. Here 
the data are smoothed in latitude with a tapered 
Gaussian having a FWHM of 10°. (A tapered 
Gaussian is a Gaussian with a quadratic function 
subtracted so that the smoothing filter and its first 
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derivative go to zero at the endpoints.) The flow is 
poleward in each hemisphere from the equator up 
to latitudes of at least 60°. The speed of the flow 
decreases as the time lag increases and the pole- 
ward flow virtually vanishes when measured with 
the 24 hour time lags. The slow down for time 
lags from 2 to 16 hour s was previously noted by 
Hathawav et al. ( 2010h . The lack of a poleward 
meridional flow for 24 hour time lags was previ- 
ously noted (and re ferred to as "anomalous") by 
Gizon et al. I ^00^. 




Fig. 2. — Meridional flow profiles for time lags 
from 1 to 24 hours. The measurements from the 
1996 data are shown with the red line while the 
measurements from the 1997 data are shown with 
the blue line. The dotted lines indicate 2a error 
limits for each year. The meridional flow is pole- 
ward up to latitudes of at least 60°. The amplitude 
of this poleward meridional flow decreases as the 
time lag increases and disappears at 24 hour time 
lags. 

Measurements with 28 hour time lags are shown 
in Figure 3. Again, the data are smoothed in lat- 
itude with a tapered Gaussian having a FWHM 
of 10°. With this longer time lag the measured 
meridional fiow is equatorward at virtually all lat- 
itudes with a peak of ~ 4 m s~^. The 2a vari- 
ations on the measurements at each latitude are 
^ 2 m s~^ making nearly all of the individual lat- 
itude measurements statistically significant. Fit- 
ting the measurements to a 4th-order polynomial 
in sinO, shown by the thick black line in Figure 3, 
gives a peak velocity of 4.6 ± 0.4 m s^^ with 2a 



errors. As such this represents the first positive 
measurement and significant (> lOcr) detection of 
the Sun's equatorward return flow. 






Latitude 

Fig. 3. — Meridional flow proflles for time lags of 
28 hours. The measurements from the 1996 data 
are shown with the red line while the measure- 
ments from the 1997 data are shown with the blue 
line. The dotted lines indicate 2a error limits for 
each year. The 4th-order fit to the combined data 
is shown with the thick, black line. 



4. DEPTH DETERMINATIONS 



Duvall I (|l98fl[ ) suggested that at different time 



lags the cells that survive are advected by flows 
with different speed (and direction) depending 
upon the depth to which those cells extend. The 
different time lags thus probe the depth depen- 



dence of the meridional flow. iHathawav I (j2012l ) 
found that supergranules of a given cell size (given 
by their wavenumber in a Fourier analysis) are ad- 
vected by the zonal flows at depths equal to their 
widths. This was determined by comparing their 
zonal velocities from the Fourier analysis with 
zonal velocities in the surface shea r layer deter- 
mined with global helioseismology bv lSchou et al. 
(Il998l) . 



Small cells have short lifetimes. Granules with 
typical diamete rs of 1 Mm have lifetimes of 



5 minutes (Title et al 



ules with typical diam eters of 
lifetime s of ^ 1 day 
iWang fc Zirin I Il989t 



Iffij 



19891 ). Supergran- 
30 Mm have 



imon & Leighton"' ^96 
Hirzbcrgc r et al. 200_ 



Cells of intermediate (mesogranular) size with 
diameters 5-10 Mm have int ermediate lifetimes 
of ~ 2 hours (iNovember et a l. 1981). As the 
time lag used in the cross-correlation increases. 
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the cells with lifetimes shorter than the time lag 
present an uncorrelated random pattern that does 
not influence the determination of the pattern 
motion. Longer time lags are associated with 
larger cells which extend deeper into the Sun. 
This assertion is supported by observations of 
the monotoni c decline in the cross-correlations 
for granules (|Title et al. I \l98^ . the monotonic 



decline in the cross-co rrelations for the superg ran- 
ule Doppler pattern ([Hathawav et al. Il2010l) .and 



the increase in lifetimes for larger supergran- 
ule cells determined fr om local helioseismology 



(jHirzherger et al. 1120081 ) 



We can determine the characteristic cell size 
at each time lag by matching the zonal veloci- 
ties from the cross-correlations with those from the 
Fourier met hod. This Fourier met hod is described 
i n det ail bv lBeck fc Schou" ( 2000l ) and Ha thawavl 
(j2012l ). Lines of mapped Doppler data at equato- 
rial latitudes are Fourier transformed in longitude 
to give spectral coefficients. Spectral coefficients 
from 6 10-day series of Doppler maps obtained at 
15-minute cadence are then Fourier transformed 
in time. The location of the peak power gives the 
zonal (longitudinal) velocity for each spatial wave- 
length or cell size. 

The line-of-sight projection for the Doppler sig- 
nal influe nces the measured zon al flow (but, as 
shown by 'Hatha way et al. ( 2010[ ) , not the merid- 
ional flow). When this projection is divided-out 
(as can only be done at the equator) the zonal 
flow as a function of wavelength from the Fourier 
analysis matches the zonal flow as a function of 
depth from global helioseismology. Figure 4 shows 
zonal velocities relative to the Carrington Rota- 
tion frame of reference for supergranules of vary- 
ing wavelength along with the zonal velocity as 
a function of dept h from global helioseismology 



(|Schou et al. Ill998f ). 



The zonal velocities for the raw data (without 
the line-of-sight projection removed) are shown in 
Figure 5. Comparing these zonal velocities with 
those obtained at the equator (also with the raw 
data) by the cross-correlation analysis gives a cor- 
respondence between time lag and characteristic 
wavelength of the cells that dominate at that time 
lag. This correspondence is given in Table 1. 




2dS0LA Rotation vs. Depth 
Supergranule Rotation vs. Wavelength 



10 20 30 40 50 60 70 80 90 100 
Wavelength or Depth (Mm) 

Fig. 4. — Equatorial zonal velocities as func- 
tions of cell wavelength and depth within the Sun. 
Zonal velocities as functions of wavelength from 
the 1996 MDI data in which the line-of-sight pro- 
jection is divided out are shown with the small 
black circles. Zonal velocities a,s functions of depth 



from global helioseismology (jSchou et al. 1 119981 ) 
are shown with the large red circles. The match 
between these measurements at virtually all wave- 
lengths and depths indicates that supergranules 
are advected by the flows at depths equal to their 
wavelengths. 



5. CONCLUSIONS 

We conclude that the Sun's poleward merid- 
ional flow is conflned to it's surface shear layer 
(from the photosphere to a depth of ^ 50 Mm). 
We draw this conclusion from the fact that the 
poleward flow weakens and then disappears as the 
cross-correlation time lag increases from 1 to 24 
hours while the measured rotation rates indicate 
that the cross-correlated patterns are dominated 
by cells with sizes, and therefor depths, that in- 
crease from ^ 7 Mm to ^ 50 Mm. A positive 



Time Lag 


Velocity 




Wavelength 


1 hour 


12± 1 m s" 


1 


7± 1 Mm 


2 hours 


24± 1 m s" 


1 


14 ± 1 Mm 


4 hours 


40 ± 1 m s" 


1 


22 ± 1 Mm 


8 hours 


49 ± 2 m s^ 


1 


26 ± 2 Mm 


16 hours 


59±2 m s- 


1 


37 ± 2 Mm 


24 hours 


71 ± 3 m s" 


1 


50 ± 5 Mm 


28 hours 


80±3 m s- 


1 


70 ± 5 Mm 



Table 1: Correspondence between time lag, equa- 
torial zonal velocity, and cell wavelength from Fig- 
ure 5. 



5 




> 2dS0LA Rotation vs. Depth 
> • Supergranule Rotation vs. Wavelength 



10 20 30 40 50 60 70 80 90 100 
Wavelength or Depth (Mm) 

Fig. 5. — Equatorial zonal velocities as func- 
tions of cell wavelength and depth within the 
Sun. Zonal velocities as functions of wavelength 
from the 1996 MDI data are shown with the 
small black circles. Zonal velocities as functions 
of depth from helioseismology (Schou et al. 199^ 
are shown with the large red circles. The zonal ve- 
locities at the equator from the different time lags 
are indicated along the vertical axis with lines indi- 
cating where they match the zonal velocities from 
the Fourier analysis. 



detection of the equatorward return flow is found 
with time lags of 28 hours corresponding to depths 
of - 70 Mm. 

The fact that the return flow signal was signif- 
icantly detected at a longer time lag in a deeper 
layer also indicates that the lack of a poleward flow 
at the 24 hour time lag is due to an actual lack of 
poleward flow rather than an insensitivity to the 
signal. 

Estimating the poleward mass flux at 30° lati- 
tude within the surface shear layer using the mea- 
sured flow speeds and the run of plasma density 
with depth, indicates that a balancing equator- 
ward mass flux at a flow speed of ~ 4 m s~^ should 
only extend another 20-30 Mm below the base of 
this surface shear layer. 

This shallow return flow is somewhat sur- 
prisi ng given earlier results from helioseismol- 
ogy- 



Braun fc Fan 



(|l998l ) used frequency shifts 
between northward and southward propagating 
waves and concluded that the meridional flow was 
poleward to a depth of at least 100 Mm 



Giles 



|2000) used time-distance helioseismology with a 
mass flux constraint and concluded that the re- 
turn flow was below a depth of about 140 Mm. 



However, iDuvall fc Hanasoge I (l2009l) have found 
that these measur ements are plagued by syst em- 
atic errors leading Gough fc Hindman I ( 2010[ ) to 
state that, with helioseismology, the meridional 
flow below a depth of 30 Mm remains uncertain. 

Recently IZhao et al. (2012A) have charac- 



terized these systematic signals, removed them 
from the data, and retrieved a meridional flow 
profile in the near surface layers that agrees 
well with measur ements obtained from magnetic 



feature motions ([Hathawav fc Rightmire 



and direct Doppler measurements (jUlrich 



■20101) 
20in[) . 



Zhao et al. I (|2012Bl ) have now reported that they 
find evidence for the meridional return flow at 
depths much shallower than the tachocline - sup- 
porting the conclusions reached here. 

A shallow return fl ow is also supported by 
Sivaraman et al. I (|2010l) who measured the merid- 



ional motions of sunspot umbrae. They compared 
the rotation rates of sunspot umbrae with rotation 
rates derived from helioseismology to find anchor- 
ing depths for umbrae of different sizes (similar to 
what was done here with supergranules) and found 
that small sunspots with unibral areas < 5/iHem 
(a /iHem is 10^^ the area of a solar hemisphere) 
are anchored at depths of ^ 90 Mm and are mov- 
ing equatorward at speeds of ^ 7 m s~^. Slightly 
deeper anchoring depths and slower equatorward 
velocities were found for larger sunspots. 

6. FURTHER IMPLICATIONS 

The shallow depth of the meridional return flow 
has implications for our understanding of both so- 
lar convection zone dynamics and solar dynamo 
theory. 

The axisymmetric meridional flow is main- 
tained by a balance between the combined ef- 
fects of three driving forces (latitudinal pressure 
gradients, Reynolds stresses resulting from corre- 
lations between radial and meridional motions in 
the cellular flows, and the Coriolis force acting 
on the axisymmetric zonal flow) and the viscous 
diffusion due to the cellular flows. Simulations 
of the effects of the Sun's rotat ion on supergran - 



ules in plane-parallel geometry ( Hathawav 19821 ) 



indicate that the Reynolds stresses produced by 
the supergranules themselves can maintain both 
the zonal shear (slower rotation at the surface 
and faster rotation below) and a meridional cir- 
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culation with poleward flow near the surface and 
equatorward flow below. The meridional Reynolds 
stresses (correlated poleward and upward motions 
in the cellular flows) are reinforced by the Corio- 
lis force on the axisymmetric zonal flow in these 
simulations. (Latitudinal pressure gradients can- 
not be maintained or included due to the periodic 
boundaries in both horizontal directions.) 

The meridional circulations produced by larger 
convection cells in rotating spherical shells extend- 
ing over several density scale heights (but thus far 
not including the surface shear layer) are usually 
highly structured in l atitude and radius a nd are 
highly variable in time (|Miesch et al. II2OO0I). More 
recent spherical shell simulations ( Miesch et al. I 
20081 ) impose a latitudinal entropy (temperature) 



gradient at the base of the convection zone to give 
rotation profiles and meridional circulations more 
in line with observations but still produce highly 
time dependent meridional flows. The source of 
this structure and variability can be attributed to 
the small number (^ 100) of convection cells that 
populate the simulated global volume. A merid- 
ional circulation driven in the surface shear layer 
by (^ 10000) supergranules (the convection cells 
that populate the surface shear layer) is far more 
likely to be less structured and variable - as is 
observed. 

The shallow return flow we find here has impli- 



individual data pairs are similar to the noise level 
from the cross-correlations. However, averaging 
the cross-correlations over as few as 24 hourly mea- 
surements could extend the results to longer time 
lags and deeper layers. Local Correlation Tracking 
using smaller patches might also yield measure- 
ments of the large scale nonaxisymmetric fiows - 
Giant Cells. 

Finally, it should be noted that, while these 
results indicate an equatorward return flow at a 
fairly shallow depth, they do not preclude the pos- 
sible existence of additional cells of poleward and 
equatorward flow at even greater depths. 

The author would like to thank NASA for 
its support of this research through grants from 
the Heliophysics Causes and Consequences of the 
Minimum of Solar Cycle 23/24 Program and the 
Living With a Star Program to NASA Marshall 
Space Flight Center. He is indebted to Ron Moore 
and Lisa (Rightmire) Upton and an anomymous 
referee whose comments greatly improved the 
manuscript and to John Beck who produced the 
temporally flltered data from the original MDI 
Doppler data. He would also like to thank the 
American taxpayers who support scientific re- 
search in general and this research in particular. 
SOHO, is a project of international cooperation 
between ESA and NASA. 



cations for flux transport dynamos (jPikpati fc Charbonneau 



Il999t iNandv fc Choudhuri I I2OO2I) - it violates a 
key assumption. These dynamo models assume 
that the return flow begins at a depth of 100 Mm 
and reaches a maximum equatorward speed of just 
1-2 m s~^ at the base of the convection zone. This 
is a critical assumption in that the equatorward 
drift of the sunspot zones is directly produced in 
these dynamo models by this slow equatorward 
flow at the base of the convection zone. 

The results presented here support the idea of 
using supergranules to probe the dynamics of the 
solar convection zone. This method can comple- 
ment the probing done with acoustic waves by the 
methods of helioseismology and provide valuable 
information on the dynamics of solar and stellar 
convection zones. The method used here, finding 
the peak in the cross-correlated Doppler velocity 
pattern for data pairs with different time lags, is 
limited to time lags less than ~ 32 hours. At 
longer time lags the correlation amplitudes from 
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